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Abstract 
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parameters using known Standard Model predictions and experimental data. 
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I. INTRODUCTION 



Mixing between a neutral meson with specific flavor and its anti-meson provides powerful 
test for the Standard Model (SM) and new physics (NP) beyond. Mixing has been observed 
l| K° — K°, B® — B®, and B® — B®. Evidence for D° mixing has been recently reported 



in 



by the BaBar [2| and Belle 0, 4] collaborations. 

Two parameters, x = AM /V and y = Ar/2r, are often used to describe the mixing 
between a meson and its anti- meson. Here T is the life-time of the meson. AM = m 2 — mi, 
Ar = T 2 — Pi with "l" and "2" indicating the CP odd and CP even states, respectively, in 
the limit of CP conservation. AM and Ar are related to the mixing matrix elements M 12 



and r 12 in the Hamiltonian by AM - zAr/2 = 2^/(M 12 - zr 12 /2)(M 1 * 2 - iY\ 2 /2). 

If a new particle has flavor changing neutral current (FCNC) interaction, a non-zero con- 
tribution to M 12 can be easily generated by exchanging this new particle in the intermediate 
state, tree or loop. The parameter T 12 must come from the absorptive part which requires 
the intermediate states be light degrees of freedom to whom the meson can decay into. 
This fact severely constrains the contributions to Ti 2 from NP. Due to this reason there is 
less theoretical work on new physics contributions to Ti 2 than that for Mi 2 . In this work, 
we study the Ti 2 parameter in the present of NP, taking SUSY R-parity violating (RPV) 
interaction as an explicit example. 

There are three types of i?-Parity violating (RPV) terms Q]: 

^fU L V L W R \ Kjk L LQi D R > ^fU-D^, (1) 

where i, j and k are the generation indices: L L , Q L , E R , D R and U R are the chiral superfields 
which transform under the SM gauge group SU(3)c x SU(2)l x U{1)y as Ll : (1,2,-1), 
E R : (1, 1, -2), Q L : (3,2, 1/3), U R : (3, 1,4/3) and D R : (3, 1, -2/3). The charge conjugated 
field ipR is defined as ip R = C{%[)r) t . We will consider each of these R-parity contributions 
to Ar 12 for meson mixing separately. In that the term proportional to \jk involves 

only leptons, it will not contribute to meson mixing, since we are not considering pairs of 
Xijk and A- jfc couplings to be non-zero at the same time. We only need to consider the last 
two terms up to one loop level. 

At the tree level, we have the following terms relevant to us by exchange s-fermions, 

C ff (\') - A ^fc A ^ , V , - / M r » f J' .J _ X ijk\*j>k> j> a j Ik ik> 
*~eff\ A ) — 2 m 2 L ^ e L u Lln u L 2m 2 U Lff1 U Lofl Ral M« R/3 



A ijk A ij'k' Jj' aij Ik jk 1 , A ijk A i' j' k -%> n i Jj' jj 
U L d R 

Kjk A 'hk'-/ a iik ,k> _ KjkKhk' a n j ik ,k> 

v hl v L a Rl^ a R r>o e LI e L a R!^ a R 



2m% L 1 L H ih> h 2m2 

d 



L 



C t J\") - ijk i,jk ' (r/ ryt* V i' J k ry _ VV' (ft jk' 

L 'eff\ A I — ^2 \ a Ra 7 a Ra a R/3! IJ- a R/3 a RaH u R/3 a Rj3! M a . 
d? 



2m 



RaJ 



'R 



A" A"* 



R 



The first two terms in C e ff(X^'^) contribute to T 12 for D° — D° mixing. Except the first 
term in C e ff(X'), all terms contribute to B% — B% mixing. 

It is clear that from the above Lagrangian at the tree level, non-zero Myi can be generated. 
Constraints have been obtained using AM for various meson mixing. However, in order to 
generate a non-zero T 12 additional loop corrections are needed from the above four fermion 
interactions. 

There are short and long distance contributions to y or T 12 . The calculations for long 
distance contributions are very difficult to handle due to our poor understanding of QCD at 
low energies. It is expected that long distance contributions become less and less important 
when energy scale becomes higher and higher, and perturbative short distance contributions 
will become the dominant one. We therefore will restrict ourselves to mesons containing 
a heavy c or b quark and to study the short distance contributions T 12 for D° — D° and 
B ls ~ B% s systems. 

For Bd y s mesons, in the SM the short distance contributions are expected to be the 
dominant ones. The prediction for Ar for B® — B® is js] 

Ar s = (0.106 ±0.032) ps" 1 . (3) 

which gives ysu — 0.078 ±0.025. The D0 experiment has measured this width difference^ 
(see also Allowing the non-zero CP violation in mixing they obtained, Ar s = (0.17 ± 
0.09 stat ± 0.03 syst ) ps^ 1 (y = 0.125 ± 0.066 sta t ± 0.022^), and in the CP conserving limit, 
Ar s = (0.12 ± 0.08 sta +° ;° 3 4s ^) ps' 1 (y = 0.088 ± 0.059 sta t t f 3 2 0syst ) . Within error bars, SM 
agrees with data. However, it is interesting to see if NP contributions can appear at the SM 
level. 

For B^ — B^ system the width difference in SM is known to be small [6] Ar^ = (26.7 ± 
0.08) x 10~ 4 ps -1 , corresponding to ysu = (2.058 ± 0.006) x 10~ 3 . There is no experimental 



data for the width difference yet. It is interesting to see whether the width difference can 
be much larger when going beyond the SM. 

n - CI 

For the D — D° mixing in the SM as was shown in [9j x and y are generated only at the 
second order in 577(3) breaking, x, y ~ sin 2 Qq x [SU (3) breaking] 2 . Most of the studies give 



x, y < 10~ 3 , although large values are not excluded [10]. 

Recently, the parameter y for D° — D° mixing has been measured. BaBar, assuming no 
CP violation in mixing, has analyzed the doubly Cabibbo suppressed (DCS) D° — > K + tt~ 
mode [2L while Belle has studied singly Cabibbo suppressed (SCS) D° — > K + K~, 7r + 7r~ 



decays [3J . From these results the authors in l_l|] have fitted the mixing parameters and get 
the following result for y with 68% and 95% probability correspondingly 

y = (6.1 ± 1.9) x 10~ 3 , y E [0.0023,0.0102]. (4) 

In this work, we find that R-parity violating contribution to the parameter y is small for 
D° — D° system, less than a few times 10~ 4 . The RPV contribution for B® — B% system 
can be larger than the SM prediction. For B° s — B° s system, the contribution to y can be as 
large as the SM contribution. 

In the following sections, we provide the detailed calculations. 



II. GENERAL EXPRESSION FOR T 12 

Before going into specific RPV model calculations, we summarize some general results 
for short distance NP contribution to Ti 2 from four quark operators generated by SM and 
NP. The calculation is straightforward. Starting from tree level four quark interactions, one 
needs to obtain the absorptive part for Fig.l. Let us take D° — D° mixing for illustration. 
For the cases considered here, we can write the AC = — 1 interaction Lagrangian as 

£ ac=-i = j2 {D M , [C x {ji)Q x + C 2 (/i)Q 2 ] + DV [C'MQz + C' 2 (fi)Q4} , (5) 
<?,<?' 

Qi = UiT^j ^TaCj , Q 2 = %Y 2 Cj , Q 3 = UiT 3 q^ qf^Ci , Q 4 = UiT^ g i r 4 c j . 

In the above we have omitted possible Lorentz indices for Tj which are contracted. The 
specific form of Tj depends on the nature of interaction generating the four quark operators. 
The notations here are that r lj2 (3,4) and r 3j 4(i )2 ) should appear on the left and right four 
quark vertices in Fig. 1, respectively. 
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FIG. 1: The one- loop Feynman diagram for meson mixing. The dashed line represents the cut for 
taking the absorptive part. 



Evaluating the diagram in Fig. 1, one obtains the following general expression for r 12 , 



1 5 



12 



^D w ,DV g (Ki^ + K 2 ^ fc ) JaO&y) (D°\ Of" \D°) 



2M T 



q,q 



a=l 



where 



(6) 



K x = (dC'iNc + (CiC 2 + C\C 2 )) , K 2 = C 2 C 2 . 



(7) 



The operators are defined as 



o\ jke = w fe r 3 7 I/ r 2 c,- UiTxY^ACi, o' l 2 



ijk£ 



Qijkl 



H fc r 3 r 2 c 3 - u e TiT 4 c h 



and the coefficients I a (x,x') are given by 

k*m 



h(x, x') 
I 2 {x,x') 



48tt 
k* 



l-2(x + x') + (x- x') 



1 + (x + x') - 2 (x - x') 



k* k* 
I 3 (x, x') = — \fx (1 + x' — x) , h(x, x') = \fx' (1 — x' + x) 



I 5 {x,x') 



k*m c 

47T 



XX' 



(8) 



where k* = (m c /2)[l — 2(x + x') + (x — x') 2 ] 1//2 with x = m 2 q /m 2 c and x' = m 2 q i / m 2 c . Replacing 
D' qq , with the SM couplings and T 3 = 7^(1 — 7s)/2 and T 4 = 7^(1 — 7s)/2, we obtain 
the formula presented in 12j, Ll3j] for SM and NP interference contribution. Note that 
when considering the contributions from the same operator, one should take D' , = D qq >, 
r"i 2 = r 34 and the expression for T 12 should be divided by 2. 



Using the above formula, one can easily work out the expressions contributing from the 
SM (taking SM operators for Qi, 2,3,4), the interference between the SM and NP (SM - NP) 
(taking Q 1>2 from SM (NP) and Q 3A from NP (SM)), and purely NP (NP - NP) (taking 
Qi, 2,3,4 from NP). New physics effects can show up in the later two cases. We will concentrate 
on these contributions. 

We comment that the fermions in the loop are not necessary to be quarks. They can be 
leptons too. If one identifies q and q' to be leptons, the correct result can be obtained by 
setting N c = 1 and Cp = in Eq. ©. One can easily generalize the above formula for 
Bj s — B® s mixing cases with appropriate replacement of quark fields and couplings. 



III. RPV CONTRIBUTIONS TO T 12 FOR D° - D° MIXING 

In this section we give expressions for contributions to with RPV interactions. 
Contributions from A' interaction to Ti2 are given by 

V2G F \' i22 \'* 12 V us V* 2 



Fl2(SM-R.PV) — 7T^ % — —X s m(Ci + C 2 )(Q) 

y y* \ 

, _ A i2k A jlk A 

12(RPV-RPV,l) ~ 19W 



y y* y y* ^2 

1 12 (RPV-RPV,l) — lno 2 2 UV/ + \Ws/J 5 l y J 



v K2rKjK'2jK*ij' m 2 c /1 

TutRPV-RPVA) - 1927rm 2 m 2 — [2™ ~ {Qs) 



where 



(Q) = {D Q \u a YPLC a u plll P LC!] \D Q ) , (Q s ) = (D°\u a P R c a UpP R c \D°) , 

The first equation in Eq. (Q is the leading order result in x s . Depending on the internal 
lepton exchanges, in the expression for r 12 (rpv-rpv,i) the indices i, j take 1 and 2 indicating 
which charged leptons are in the loop. In principle, one can also have an electron and a 
tauon in the loop. However, the tauon mass is close to the D meson mass, the contribution 
is suppressed by phase space. We will neglect this contribution. In the expression for 
I" '12 (rpv- rpv, q)i hf take 1 and 2 indicating which of the down quarks are in the loop. 

Fi2(sm-rpv) comes from SM interaction with the second term, Ti 2 (rpv-rpv,i) comes from 
the first term, and Ti 2 (Rpv-RPV,q) comes from the second term, in £ e //(A'), respectively. 
Note that the SM-RPV contribution is proportional to the internal quark masses and the 
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dominant one comes from ss in the loop. This is due to the chiral structure of Tj which 
allow only 0$ kl to contribute and therefore proportional to the function I$(x,x'). If x or x' 
takes the down quark mass I$(x,x') is negligibly small. 

In obtaining the expression for ^v2(sm-rpv)i we have used the SM AC = — 1 Lagrangian, 

4G F 

C SM = -^j=-V* c V q , u [C x {m c )Q x + C 2 (m c )Q 2 ] , (10) 

with Ti and T 2 in Eq.© to be 7^(1 — 75) /2 and 7^(1 — 7s)/2, respectively. 

The contributions from A" interaction come from the first term in £ e //(A") and are given 

by 

r 1KSM - RPV) = -^%^%^x s m^[(2C 1 + C 2 ){Q>) - C 2 (Q'}} , 

d R 

T^rpv-rpv) - 1927Tmlml/ — ^W)J ' ^) 



4 4 



where 



(Q') = {D Q \u a rPLC a u pl ^P R c p \D°) , (Q') = {D \u a ^P L c p u pl ^P R c a \D G ) . 



In first equation of Eq. ffTTj) . as in the first equation of Eq. (J9J), we only kept the leading 
order in x s . The SM-RPV contribution is dominated by ss pair in the loop for the same 
reason as that for the A' case for SM-RPV contribution explained earlier. 

Here we should mention that recently in Ref. [3] the authors have considered RPV 
with slepton and squark exchanges for SM-NP contributions. Our predictions in the first 
equations in Eqs. ([9]) and (llip . for the same measurable, do not agree with their Eqs. (16) 
and (24), respectively. 

IV. RPV CONTRIBUTIONS TO T 12 FOR - B° d s MIXING 

In this case all terms except the first term in £ e //(A') contribute to Ti 2 . 

A. The A' contribution 

The expressions for Ti 2 from various contributions are given by 

^(sm-rpv) = ^ FmlX ^ k {(2dK) - C 2 {m b ))(Q') + (2C 2 (m 6 ) - C 1 (m 6 ))(Q')} , 
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_ m b } A j3i' A j'ki' A j'3i A jki n ^ , 

1 12(RPV-RPV,v) — 77^ \ 2 2 W + \Vs/J 

V d R d R 

+ ((g) + (q.» - 2 ^4^ ^# f(Q') - l(g' 

777. 7 777. ~ ./ 777. ~ / 777- r- \ z 

1 12(RPV-RPV,l) — 77^ 2 2 UVs/ + \V/J > U^J 



m 



\ 1 \i* \/ \/* , -1 

A ijk A ij>3 A i'j'k A i'j3 I 1 



?12 { RPV-RPV, U) ~ 192nmB m% [^(Q) ~ (Q 



Fl2(RPV-RPV,d) ~ 1927r b mB m \ m ^_, { 18X 'ijj' A ^3 A 'i'3kK*jj'(Q') 

+ {K3j'K*kjK'3jK*kj' + KjkK*j'3K'j'kK*j3) ~~ 

where j,j' take the values 1 and 2. 

Accfc = Kg fe ^cbA- 2A ,A-23 , A cu fc = V* qk V u bX' ilk X'^ 3 , 



A uc fc — V* qk V cb X' i2k X'* l3 , A UMfc — V^V^A^A^ . 

and 

(Q') = (B° qk \q k a ^P L b a q h ^P R bp\B° qk ) , 
(QO = (Sj^^/J qfaP R b a \B qk ) . 

Here for and B s systems, k takes 1 and 2, respectively. 

The five different contributions to Ti 2 listed above come from the second, first and sixth, 
fourth and fifth, second and third terms in C e ff(X'), respectively. 

B. The A" Contribution 

In this case we have 

V2G F x c VT^4r c mlV; qk V cb K kl X' 2 % 

1 12 (SM-RPV) — g 2 



a R 

{(20^) + C 2 (m b ))(Q') - C 2 (m b )(Q>)} , 



m 



2 



\ ii \ri* \ ii \ii* , q 
A ijk A i'j3 A i'j'k A ij'3 I 6 



r i2{ RPV-RPV,u) ~ lg o ^Q) ) > 

d R d R 

2 \il \ll* \tf \ II* , Q 

m b A j'ki A j'3i' A jki' A j3i ( 6 



1 12(RPV-RPV,d) ~ 77 2 2 o\V/ 



w 



R 
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where i, %' take the values 1 and 2. The first two terms are due to the first term in C e ff{\"), 
and the last term is due to the second term in C e ff{\")- 

The SM-RPV interference is dominated by cc exchange in the loop for the same reasons 
as that for the ss dominance for SM-RPV D° — D° mixing. 



V. NUMERICAL ANALYSIS 



In this section, we carry out numerical analysis for RPV contributions to mixing param- 
eter y for D° — D° and B® — B^ s systems. 

In general RPV contribution to Ti2 has CP violating phases associated with the new 
couplings. The relation of Ti2 and y is not trivial. If CP violating effects can be neglected 
which is the case for the SM, they have a simple relation 

Ti2 ,.. t \ 

V=jr- (14) 

In our numerical analysis, we assume CP conservation for easy comparison with data and 
other constraints obtained in the literature. 

To compare with data, one needs to evaluate various hadronic matrix elements in the 
expressions for T^. We write them in the following form 

(Q) = \f P m P B Q , (Q s ) = -^f P m 2 P B Qs , 

(Q') = -\f P m P B QI , <Q'> = --f P m P BQ, , (15) 

where Bq factors are the so called bag parameters [ljj] . This way of parameterizing the ma- 
trix elements was inspired by vacuum saturation approximation. In the vacuum saturation 
approximation, they are all equal to one, which we will use in our estimate. 

In the Table I we list the parameters and coefficients appearing in the equations above. 



The input CKM elements are 16] 



A = \V US \ = 0.2248, AX = \V cb \ = 41.5 x KT 3 . (16) 

The charm quark mass also comes into the calculations. In our numerical analysis we identify 
m c and nib with pole masses. Numerically we use jl4] m c ^ i e j m\,^ i e = y/x~ c = 0.29, which is 
based on the mass difference m CiPO ; e — nibble = 3.4 GeV and mb ;PO ie = 4.8 GeV. 

To give some understanding of RPV contributions, in the following analysis we take the 
central values for the input parameters. 



H = 1.3 GeV [i = 4.8 GeV Masses (GeV) Decay cons. (GeV) Widths (GeV) 
[6,14] [6,14] fl] [15J fl] 



x s = 0.006 x c = 0.0841 M D = 1.8645 
Ci = -0.411 Ci = -0.272 M Bd = 5.279 
C 2 = 1.208 C 2 = 1.120 M Bs = 5.368 



f D = 0.201 r D = 1.6 x 10" 1 
/ B = 0.216 r Bd = 4.27 x 10- 
f Bs = 0.260 T Bs = 4.46 x 10- 



TABLE I: The central values of input parameters and coefficients. 



A. D° — D° mixing 



Taking all A' and A" to be real, and inserting known values for the parameters involved, 
for A' contributions, we have 

(lOOGeV) 2 



V(sm-rpv) = 0.0037 x A- 22 A 



i* 

12" 



2 



V(rpv-rpv,i) - 0.3298 x \ 2 kK*ik^j2k'Kik' 



(lOOGeV) 4 
m k m k 



(17) 



- 9893 x \' \'* V A'* ( 100Ge ^) 4 
U(RPv-RPv,q) — u.y»yd x \2j' A nj A i'2j A i'ij'~^2 



The A" contributions are given by 

V(SM-RPV) : 
U(RPV-RPV) 



-0.0077 x A^ 2 A 2 '* 2 



(lOOGeVf 



d R 



1.3191 x AyjA 2 ^Ay V A 2 */j 



{WOGeVy 



(18) 



17, HQ- 



There are constraints on the RPV parameters from various other processes 
Taking these constraints into account, we list in Table II the corresponding values for the 
mixing parameter y. 

For the contribution due to the A' terms, using the constraint |A- 2fc Aj* lfc | ^ 5.28 x 10~ 6 , 
obtained from K + — > tt + i'u H] and taking the same bound for \X'i2j^ij\j=i,2 with the 
assumption that there is no accidental cancellation, we find V(rpv-rpv,i) an d V(RPv-RPV,q) 
to be tiny < 10 -11 . For the interference term due to the SM with A' term, we have not 
found direct constraint on the appropriate combination of the A' terms. We therefore use 
individual constraints from (Ai 22 Ai* 12 , A 222 A 2 * 12 , A^A^) = (0.0009,0.0124,0.0572). This 
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RPV parameters Bounds [Processes] 



Estimate 



\ A i22 A il2\ 


0.07 (see text) 


V(SM- 


-RPV) 


~ 2.6 x 10" 


i 


\ A i2k A jlk\ 


5.28 x 10" 6 [K + -► tt h 


~W] V(RPV- 


RPV,l) 


~ 9.2 x 10" 


-12 


\ A i2j A i*lj\j=h2 


5.28 x 10" 6 [K + -> vr^ 


'vv\ y(R PV - 


RPV,q) 


~ 2.5 x 10 


-11 


l A 132 A 232l 


3.1 x 10" 3 [£>£>] 


V{SM- 


-RPV) 


~ 2.4 x 10" 


5 


l A 132 A 232l 


3.1 x 10" 3 [£>£>] 


U(RPV 


-RPV) 


~ 1.3 x 10 


-5 



TABLE II: The bounds on parameters from 



13, 



18 



19, 



2l| and corresponding values for y. 



leads to D(sm-rpv) — (3-3 x 1CT 6 , 4.6 x 1CT 5 , 2.1 x 10~ 4 ). To see the largest possible value for 
y, we sum these three with the same sign to obtain an upper bound V(sm-rpv) = 2.6 x 1CT 4 . 
This contradicts with the result obtained in Ref. [12I . where y can be as large as ~ —3.7%. 



As for the contributions from A", from the constraint [A'^A^I £ 3.1 x 10 3 



2l| , 



A'/ 12 , A'/ 2 i ~ 2 x 10 9 and \" 31 < 10 4 [18], we find that A" contributions are small, less than 
TO" 4 . 

We conclude that R-parity violating contributions to y for D° — D° mixing are small. 
B. B° d - B% mixing 

Here we present the numerical results for the — B d mixing. For A' contributions, we 
have 

00, = x (WOGeV) 2 
U(sm-rpv) = —224.5 x X qlql ^ ; 



V(rpv-rpv,u) = 55.1 x (lOOGeV) 



4 J A j3i' A j'li' A j'3i A jli 

2 2 

m%, m% 

d R d R 



A jil A j'i3 A j'i'l A ji'3 , r } A j3i' A j'li' A j'il A ji3 

2 2 2 2 ( ' 

U{RPV-RPV,l) — OO.l X A j ii A ji i 3 A j'i'l A ji l 3 2 2 ' l ly J 

V(RPV-RPV,u) — i-OO.Z X A iji A ij^ A i/jii A i/j 3 2 2 ' 

TYl TYl ! 

y{RPv-RPv,d) = 165.2 x [K3j'K*ijK'3jK*ij' + KjiK*j'3K'j'iK*j3 (20) 
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RPV parameters Bounds [Processes] 



Estimate 





1.2 X 10~ 5 


[BB] 


U(SM- 


-RPV)! 


~ 1.1 x 10~ 4 


^21^23 


5.0 x 


[SB] 


V(SM- 


-RPV) 2 


~ 1.0 x 10 -4 




1.1 x 10 6 [B u 


■xr — 1 


U(RPV- 


-RPV,u) 


~ 2.7 x 10~ 4 




l.l xlO j [B u 




U(RPV- 


-RPV,l) 


- 0.67 x 10~ 4 


A ijl A ij'3\ ' \ A i'j'l A i'j3\ 


6.4 x 10~ 7 


[BB] 


U(RPV- 


-RPV,u) 


~ 1.1 x 10~ 4 


x' v 1 1 \> \' 1 

A ijj' A il3\ \ A ijj' A i3l\ 


1.6 x 10~ 6 


[BB] 


V(RPV- 


-RPV,d) 


~ 7.2 x 10~ 3 


l A 212 A 232l 


6 x 10~ 5 [B 




V(SM- 


_ RPV) ~ 2.8 x 10- 5 


l A i21 A j23l 


6 x 10~ 5 [B 


-> 07r] 


V{RPV- 


-RPV,u) 


~ 0.8 x 10~ 6 


l A j'12 A j32l 


6 x 10~ 5 [B 


-> <H 


V{RPV- 


-RPV,d) 


~ 3.2 x 10" 6 



TABLE III: The bounds on parameters from 

y- 



17 



18 



19 



22 



24 



251 ] and corresponding values for 



00 \' \'* \' \'* 

Z ° A ijj' A il3 A i'31 A i'jj' 



(IQOGeVy 



For A" contributions, we have 



V(sm-rpv) — 0.46 x A212A232 



(lOOGeVf 



990 * x A" A"* A" A"* ( 100GW ) 4 
D(RPV-RPV,u) - ^U.d X \21 A i'23 A i'21 A i23 ™2 



m h m h 



881 q v x" x" V'* (lOOGeV) 

V(RPV-RPV,d) — ooi.O X A j , 12 A j , 32 A jl2 A j32 



m 2 ,,m 2 



(21) 



We list various constraints on relevant RPV parameters and corresponding values for y in 
Table III. 

F° r V(sm-rpv), we keep only the two terms proportional to A cc i and X uc i since the other 
two terms are proportional to V u ^. We obtain, A g/9 i = v4A 2 A- 21 x [A^* 3 — AA- 23 ]. In Table III, 
ViSM-RPV)! and U(sm-rp) 2 . indicate contributions from the first and the second term in A g/g i. 
Using constraints from [29J, we have |A- 21 A£ 3 | < 1.2 x 10~ 5 and |A- 21 A- 23 | ^ 5.0 x 10" 5 , each 
gives y ~ 1 x 10~ 4 . This value is much less than the SM prediction. 

For y(Rpv-RPV,v{i)), using the constraints lA^A-,^], |A-yA-, 3j | < 1.1 x 10 -3 , from [19] we 
find that the corresponding upper bounds: y(RPv-RPV,u) — 2.7 x 10~ 4 and V(rpv-rpv,i) — 
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6.7 x 1(T 5 . 

As for the contribution U(rpv-rpv,u), there are four terms with take values 1 or 2. 
Taking the explicit constraints from Ref. [22j, |A- U A- 13 | < 8.0 x 10~ 4 , |A. n Ay < 2.5 x 10~ 3 , 
|A^ 21 A^ 13 | < 1.2 x 10~ 5 , and |A^ 21 A^ 23 | < 5.0 x 10~ 5 , we find that the dominant contribution is 
from the case j = j' = 1 which gives the upper bound y(Rpv-RPV,u) — 1.1 x 10~ 4 . In the same 
way for the contribution y(Rpv-RPV,d), taking the constraints from [22j, |24j, the dominant 
part is from A^A^A^A^-/ term with j = j' = 1. We find the value for y(Rpv-ppv,d) can 
be as large as 7 x 10~ 3 . This is about three times larg er than the SM contribution. 



Contributions from A" are also constrained. Ref. |25[ considers the decay mode B~ — > <f)7r~ 
and drives the upper bound A" 21 A" 32 < 6 x 10~ 5 . Taking the same bound for |A" 21 A" 32 |j =2 
and |A" 21 A" 32 |j = i j2 under the assumption that there is no accidental cancellation, y(sM-RPV) 
and y(Rpv-RPV,u{d)) are constrained to be small as can be seen from Table III. 

We conclude that if there is no accidental cancellations, for B® — B® mixing, R-parity 
contribution to y can be as large as 7 x 10~ 3 which is about three times larger than the 
SM prediction. This value is still difficult to be measured experimentally. However, if there 
is accidental cancellation, y could be bigger. Careful measurement of y for B® — B® can 
provide valuable information about new physics beyond the SM. 



C. B»-BU mixing 

For A' contributions, we have 



x (lOOGeV) 2 
V(sm-rpv) = -316.6 x \ q , q2 2 ' 



y(RPV—RPV,u) 

77.7 x (lOOGeV) 



4 J A j3i' A j'2i' A j'3i A j2i 
2 2 

m%, m% 

d R d R 



A ji2 A j'i3 A j'i'2 A ji'3 . A j3i' A j'2i' A j'i2 A ji3 
' 2 2 ' 2 2 I ' 

m 4 ™% m di 

— 77 7 x A' A'* A' A'* ( 100GeV ^ (22) 

y(RPV-RPVJ) — It-lX A ji2 A j'i3 A j'i'2 A ji'3 2 ^2 ' \ ZZ > 



, v x' x- x' x- ( 100GeV ) A 

y(RPV-RPV,u) — ±o6.L X \j2 A ij>3 A i>j'2 A i>j3 ^2 ™2 ' 



my mi, 



y(RPv-RPv,d) — 233.1 x (A- 3j /A- 2j A-/ 3j A-T 2j -/ + A- J - 2 A i */ 3 Aj/ J / 2 A i T ? - 3 
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98 X' X'* X' X'* A (lOOGW) 4 

For A" contributions, we have 



9 Q v X" X- ( 100GW )' 
U(SM-RPV) — —^■v X A 221 A 231 2" 



U(rpv-rpv,u) - 310.8 x A" 12 A "* 13 A ", 12 A "*3 - 2 2 - , (23) 



"a "a 



-1943xA" X"* X" (lOOGeK) 4 

y{RPV-RPV,d) — i24<j X A i'21 A i'31 A i21 A i31" 



RPV parameters 



Bounds [Processes] 



Estimate 



Our bounds on RPV 



I X' X'* I 
\ A i23 A i22\ 

I \/ \/* I I \ /* \/ I 

l A i3j' A i'2i'l> \ A ij2 A i'j3\ 

l A ij2 / Vj3l 
l A ij2 A ij''3li,iY3 

I A i23 \'32 '+1 

\//* I 
l A 221 A 23ll 

l A il2 A il3U^3 
l A j'21 A i'3ll 



8.2 x KT 3 [B S B S ] y(SM-RPV) - 0.11 7.4 x KT 3 (7.8 x 10~ 4 ) 

1.5 x 10- 3 [B A>i>] y^py-W,.) - 7 x 10~ 4 

1.5 x 10- 3 [B - X s ^] j/^v-jiPv.o - 1-7 x lO" 4 

5.16 x lO" 2 [J3 a J a ] y(RPV-RPV,u) ^ 0.06 2.0 x 10~ 2 (6.6 x 10~ 3 ) 

see text V(RPV-RPV,d),3 - 0-26 1.5 x 10~ 5 (1.6 x 10~ 6 ) 

1.01 x 10~ 2 [B -» Kit] V(sm-rpv) - 2.9 x 10~ 2 3.3 x 10~ 2 (3.5 x 10~ 3 ) 

see text V(rpv-RPV,u), see text 1.77 x 10 -2 (5.7 x 10~ 3 ) 

1.2 x 10- 3 [B+ -> ir+K°] y iRP v-RPv,d) - 1-8 x 10~ 3 



TABLE IV: Upper limits on parameters from 



17 



18 



22, 



23, 



24, 



26J and corresponding values 



for y. The numbers in the brackets correspond to the case, when central values for ysM and yExp. 
are used to put the constraints. For each number see the text for the explanation. 



HE? 



We list the constraints on the RPV parameters from , LL9|] and the corresponding values 
for the mixing parameter y in Table IV. 

There are several terms contributing to y from A'. For y(SM-RPV) case we again drop 
terms proportional to V u b, and have, X q/q 2 = v4A 2 A- 22 x [-M-a + ^-Mis]- We are using constraints 



from Ref. 



23| we have \X' m X'^ 2 \ < 8.2 x KT 3 and from Ref. |24j |A< 13 A£ 2 | < 2.48 x KT 3 . 



The first term dominates and gives V(sm-rpv) — 0.1, which is of order of SM prediction 
y SM ~ 0.078 and may have measurable effect. 

For nRPV ^„ we have three contribution, For first and second contentions using 
the following conditions on RPV parameters |A'- 3i /A'* 2 j/|, |A-* 2 A-/ j3 | < 1.5 x 10 3 [17], we get 
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U(rpv-rpv,u) — 1-7 x 10~ 4 . For the last term, we obtain V(rpv-rpvm) — 3.5 x 1CT 4 . If we 
simply add them together we will get D(rpv-rpv,v) — 7 x 1CT 4 . 

For y{Rpv-RPV,i)i the situation is the same as the second term of {RPV — RPV, y) case. 



In the case for V(rpv-rpv,u), if one uses the individual constraints from 19j, |23j, |24] 
(A^^AtaASs^A&.AtoASa) = (1.63 x 10- 3 ,8.2 x 10~ 3 ,5.16 x 10- 2 ,2.48 x 1(T 3 ) can 
get for each contribution U(rpv-rpv,u) — (6.2 x 10~ 4 , 1.6 x 10 -2 , 6.0 x 1CT 2 ). If we keep the 
dominant interference term we will get y ~ 0.06. 

For y(Rpv-RPV,d) case we have three contributions. The contribution of the first term 
into y is small, about 1.6 x 10~ 3 . The dominant contributions here are coming from squares 
of l A i3i A Sil - L29 x 10-3 \ B ~ -" K ~^o] fl and |A^ 32 A£ 2 | ~ 2.3 x lO" 3 [B° -> MM] 
[l9| . The second term is just the same as in y(Rpv-RPV,u) case, considered above. So 
here we have y(Rpv-RPV,i) — 0.06. The last term is enhanced with the large coefficient. 
Here the dominant contributions are coming from |A^ 22 A-2 3 A^3 2 A-T 22 | ~ 18.9 x 10 -6 and 



|A: i2 A£ 3 A;, 32 A^ 12 | ~ 20.6 x 10- 



22 



23 



241 ]. If we simply add them together, their 



contribution will be y(Rpv-RPV,d) — 0.26. So, here we can conclude, that in y(Rpv-RPV,d) 
case one also can expect large effects for y. 

For A" contribution to ytSM-RPV), we have |A 221 A 231 | ^ 1.01 x 10~ 2 from B — > Kti |24j, 
which gives y(sM-RPV) — 2.9 x 10~ 2 . 

For the y(Rpv-RPV,u) case no direct constraint onJAf 12 A-' 13 |j = i i2 exists. However, if one 
assumes that |A? 12 A£l,| i=lj2 w |Af 12 A^ 3 | <1.2 x 10" 3 [2fl|, y {RP v-RPV,u) ^ 4.5 x 10~ 4 . 

For the last case of Eq. ( 123]) from [26( we have |A", 21 A", 31 | < 1.2 x 10 -3 , which gives 
y(RPV-RPV,d) ~ 1-8 x 10~ 3 . 

We note that present constraints on the RPV parameters still allow large y(sM-RPV), 
y(RPV-RPV,u) and y(Rpv-RPV,u) from A' interaction. One can turn the argument around 
to constrain the relevant RPV parameters by requiring that the new contributions do not 
exceed the allowed range for the difference of SM prediction and experimental values. We 
have carried out an analysis, taking the one sigma range values ysM = [0.054,0.101] and 
VExp. = [0.022,0.151], and assumed constructive contributions between the SM and new 
contributions, to obtain the bounds for each individual terms. Similar analysis has been 
performed for A" cases. The bounds are listed in Table IV in the last column. These bounds 
are new ones. 
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VI. CONCLUSION 



In this paper we have explored the influence of SUSY R-parity violation contributions 
for the lifetime difference y on the D° — D° and B^ s — B® s systems. We have obtained 
general expressions for new physics contributions to y from effective four fermion operators 
including SM-NP interference and pure new physcis contributions. We find that in general 
R-parity violating contribution to D° — D° mixing, and B® — B® to be small. There may be 
sizable contribution to B° s — B® mixing. We also obtain some interesting bounds on R-parity 
violating parameters using known Standard Model predictions and experimental data. 
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